Spin-Hall conductivity of a disordered 2D electron gas with Dresselhaus spin-orbit 
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The spin-Hall conductivity of a disordered 2D electron gas has been investigated for a general 
spin-orbit interaction. We have found that in the diffusive regime of electron transport, the dc 
spin-Hall conductivity of a homogeneous system is zero due to impurity scattering when the spin- 
orbit coupling contains only the Rashba interaction, in agreement with existing results. However, 
when the Dresselhaus interaction is taken into account, the spin-Hall current is not zero. We also 
considered the spin-Hall currents induced by an inhomogeneous electric field. It is shown that a 
time dependent electric charge induces a vortex of spin-Hall currents. 
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PACS numbers: 72.25.Dc, 71.70.Ej 

Spintronics is a fast developing area using the electron 
spin degrees of freedom in electronic devices^* 2-3,4 . One 
of the most challenging goals of spintronics is to find a 
method to manipulate spins by electric fields. The spin- 
orbit interaction (SOI), which couples the electron mo- 
mentum and spin, can serve as a spin-charge mediator. 
There have been several suggestions to use the SOI in 
semiconductor quantum wells (QW) to create the elec- 
tron and hole spin currents and to accumulate the spin 
polarization by applying an electric field parallei£*£iL& or 
perpendicular^iS to the QW. The spin current induced 
by the parallel electric field and flowing perpendicular to 
it has been named the spin-Hall effect (see alsoAi) . Since 
the prediction of this effect by Murakami et. alA and 
Sinova et. al&, there have been much discussions con- 
cerning the effect of nonmagnetic impurity scattering on 
the spin-Hall conductivity in systems with Rashba spin- 
orbit coupling. Some groups predicted that the impurity 
scattering should suppress the spin-Hall effect induced 
by a homogeneous and static electric fieldi2ii 3 *± 4 *i& ; even 
if the mean scattering time r is much longer than 1/A, 
where A is the spin-orbit splitting of the electron energy 
(we set h=l). This result was confirmed by an analysis 
of the sum rules in Refii^. Yet some other groups came 
to different conclusionsiii&iS. 

In the present we use the diffusion approximation to 
derive an expression of the spin-Hall conductivity for 
a general SOI including both Rashba and Dresselhaus 
terms. For pure Rashba SOI, as well as for linear Dres- 
selhaus interaction, we found that the dc spin-Hall con- 
ductivity of the homogeneous system becomes zero even 
for a weak disorder scattering, confirming thus the re- 
sults of Ref. [T^ITMi3h1 Hi . On the other hand, when 
the cubic terms of Dresselhaus SOI is included, a finite 
spin current is produced. In order to study the effect 
of a spatially inhomogeneous electric field, our analysis 
keeps finite frequency O and wavenumber Q of the elec- 
tric field. We found that for Q,<^DQ 2 , where D is the 
electron diffusion constant, the flow of the spin-Hall cur- 
rents is dominated by the screening effects. Similar to 
formation of an electron screening cloud around an ex- 



ternal charge, the spin-Hall currents form a vortex. 

We consider a typical III-V semiconductor QW with 
only the lowest subband occupied. The spin-orbit cou- 
pling of conduction electrons has the form 



H, 



(1) 



where <r=(a x , (J y ,a z ) is the Pauli matrix vector, and 
a function of the two-dimensional wave- vector k. In gen- 
eral, hk contains both the Dresselhaus and the Rashba 
terms. The former exists also in bulk crystals^, while 
the latter appears only in asymmetric QWs^i. For a 
QW grown along the [001] direction, which is set as the 
z axis, the Dresselhaus SOI is given by 2 ^ 



h% = 
K = 



Pk x (k 2 y -a 2 ), 
-(3k y (k 2 x -a 2 ), 



(2) 



where the parameter a 2 is the average of the operator 
-(d/dz) 2 with respect to the lowest subband wave func- 
tion. The Dresselhaus SOI in J2J contains terms both 
linear and cubic in k. Usually, in heavily doped QWs, 
for electrons at the Fermi energy both terms are of the 
same order of magnitude 2 ^. The Rashba interaction has 
the form 21 



h x 



ak y 



K 



-ak x 



(3) 



Let us apply an electric field along the x axis, and 
express it as the gradient of a scalar electric potential 
2£=-VV. This gauge is more convenient for studying 
the case of finite wave-numbers Q in the Fourier ex- 
pansion of E. The one-particle spin-current operator is 
jj=(CT*u J +t; J <7 t )/4, where the particle velocity is 



_^ + A(ft fc cr) 

to* dk % 



(4) 



This definition has to be used with cautious, since the 
spin current is not conserving in systems with SOI, as 
discussed in Ref. 0. We are interested in calculating 
the spin current polarized in z direction and flowing in y 
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direction. Since h k in J3J and has no z components 
the spin-current operator is Jy=cr z k y /(2m*). We will 
calculate the corresponding spin-Hall current within the 
standard linear response theorjiS 5 and denote it as J. So, 
the initial expression for J is 

k.k' 71 ^ 

xJ*G r (k' + ,k + ,u + n)])V(n,Q), (5) 

where k±=k±Q/2, and n F {oj) is the Fermi distribution 
function. In JSJ the trace runs through the spin variables, 
and the angular brackets denote the average over the ran- 
dom distribution of impurities. The terms containing the 
products of the form G a G a and G r G r are neglected since 
their contribution to the spin-Hall current is small— For 
simplicity we assume that in the vicinity of the Fermi en- 
ergy Ep , the amplitude of impurity elastic scattering is 
isotropic and momentum- independent. In the quasiclas- 
sical approximation, when E f t^>1, the average of the 
product of the retarded and advanced Green functions 
G r and G a can be calculated perturbatively. If we ignore 
weak localization effects, the perturbation expansion of 
(JSJ) consists of the so called ladder diagrams 25-26 . For 
small and Q these diagrams describe the particle and 
spin diffusion processes. The spin diffusion also includes 
the D'yakonov-Perel spin relaxation— Therefore, the 
spin-Hall current JSJ) is determined by the combination 
of spin and particle diffusion propagators. 

To calculate and to combine these propagators forar- 
bitrary h k , we will follow the formalism of Ref. |28l29| . 
In JSJ the spin-current vertex J z is coupled to the spin in- 
dependent potential V. Such a spin-charge coupling has 
two channels. In the first channel, J z and V are coupled 
via the spin-independent particle diffusion propagator. 
This contribution to the spin-Hall current is denoted as 
J\. For Q<^1/t and v f Q-^1/t, where v F is the Fermi 
velocity, from JSJ we obtain 

Ji = i^D(n,Q)v(n,Q), (6) 

where D(Cl,Q)=[r(-i£l + DQ 2 )}^ 1 is the particle diffu- 
sion propagator—. The vertex \I> is 

* = £ Tr[G r (k + , E F + n)G a (k_,E F )J*} , (7) 

fc 

where G r,a (k,E) are the Green functions averaged over 
random impurity positions. 

The second coupling channel is more complicated. The 
spin current couples first to the spin diffusion-relaxation 
propagator, which couples to V via the mixing of charge 
and spin diffusion processes. The mixing of these diffu- 
sion processes was pointed out explicitly by Burkov at. 
ai— . The spin-Hall current due to this channel is de- 
noted as J2, and is obtained as 

J 2 = i^ l D l i(n,Q)M j D{Sl,Q)V(Q,Q), (8) 



with the vertices 

¥ =Y / Tr[G r (k + ,E F + n)<T l G a (k_,E F )J*}. (9) 
fc 

In JHJ the superscripts I and j are summed over x, y and 
z. The spin diffusion-relaxation propagator Z) lJ (fi, Q) 
describes diffusion and relaxation of a spin density 
packet. Therefore, this propagator satisfies the spin diffu- 
sion equation for spins polarized in the j-direction when 
a source creates spins polarized in the i-direction. M 3 is 
the spin-charge mixing, defined as 

Mi - ^-X Tr[Gr{k+ > EF+n)x 

u fc 

xG a (fc_,£ F )^"], (10) 

where N =m* j '(2ir) is the 2D density of states. M° 
makes the diffusion of spins polarized in j-direction de- 
pendent on the charge density distributioni&it This 
spin-charge coupling is weak and is proportional to the 
small parameter h k /E F . Therefore, in (JSJ) we keep only 
the terms linear in M° . It should be noticed^ that J2 is 
closely related to the electric field induced accumulation 
of the in-plane polarized spin density S . For example, it 
can be shown that J 2 =^ 1 S l /2tttN . 

After averaging over the impurity positions, the re- 
tarded and advanced Green functions are obtained as 

G r {k,E) = [G a (k,E)}^ = (E-Ek-hk-a+iT)- 1 , (11) 

where r=l/(2r) and E k —k 2 /(2m*). For the case of 
short-range impurities and the constant density of states 
near Ep, the scattering rate T is independent of momen- 
tum— Using Qllfl. for small f2 and Q, one gets from |J7J, 
© and EU 



¥ = -ttN e lmz vyh^Zh , 
M j = ^Q m (V£X)/i|Z fc , (12) 

where Z k =(T 2 +h k )~ l and n k =h k /h k . The over-line in 
(|12fl denotes the average over directions of fc which has 
the magnitude k=k F . In (|12l) e lmz is the antisymmetric 
tensor with e xyz —l, and all doubly repeated superscripts 
should be summed over x, y and z. 

D y '(f2, Q) satisfies the spin diffusion equationi^^S. For 
Qvp<^h kp we can neglect in this equation the diffusion 
and spin precession terms which are proportional to the 
gradient of the spin propagator. We then have 

-iVLD m3 \n, Q) = 2T5 mj - T ml D l3 (Q, Q) , (13) 

where T rnl is the spin relaxation matrix element. At 
low frequency the relaxation term dominates and so 
D m3 (Q, Q) is simply given by the inverse of T ml , and 



T ml = 2T[S ml h 2 -h%h l k }Z k . (14) 
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This equation differs by a factor T 2 Zk from the standard 
definition of the spin relaxation matrix, for example, in 
Ref. |23- This factor is not unity because we consider 
the situation that the spin splitting A=2/ifc can be com- 
parable to the electron elastic scattering rate 2Y. 

Let us first consider the case of Rashba SOI ©■ We 
then set Q y =0 and E=-iQ x V to calculate $> y and M y 
from l|12|) . In this case both the spin relaxation matrix 
and the spin diffusion-relaxation propagator are diagonal. 
Substituting the so calculated H y , M y and D yy into 
© and |[SJ, the currents J\ and J 2 are obtained as 



Ji =-J 2 = E 



n 



8tt 4F 2 + A 2 n + iDQ 2 



(15) 



where A=2akF ■ Hence, the total current J1+J2 vanishes 
even for small impurity scattering rate T<cA, in agree- 
ment with the existing resultsi2iiiiiiSii&. We should 
mention that in deriving this result for Cl^X yy , in the 
denominator (-i£l+T yy ) of the spin diffusion-relaxation 
propagator the frequency term has been removed. If we 
retain fi, J\ and J2 will cancel each other not exactly but 
the accuracy^ is up to Q,/T yy . As was pointed out by 
Mishchenko et. ali&, near the sample boundaries J2 can 
also differ from J\ because of the rapid spatial variation 
of the spin diffusion propagator. We have ignored this 
effect by neglecting the gradient terms in the diffusion 
equation i|13|) . If necessary, in our approach we can con- 
sider the boundary problem by substituting into © the 
complete solution D m i (fl,Q) of the spin diffusion equa- 
tion^ Our main goal is, however, to show that the spin 
current is not zero in the bulk of the sample when the 
Dresselhaus SOI is taken into account. In this case the 
total spin accumulation near the sample edge will be de- 
termined by a direct inflow of the spin polarization from 
the bulk. 

Let us assume that the SOI contains only the Dres- 
selhaus interaction J5J, which has terms both linear and 
cubic in k. When the cubic interaction is ignored, there 
is no spin-Hall effect because the linear Dresselhaus SOI 
can be obtained from the Rashba SOI via a unitary trans- 
formation of the spin operators^. For the complete Dres- 
selhaus interaction ©, following ©, © and lfT2|) - (fil|l . 
the calculation of the spin-Hall current is straightforward. 
We obtain the total spin current J= Ji+ J2 as 



J = Ea, 



H 



n + iDQ 2 



(16) 



where ct s h is the DC spin-Hall conductivity at Q— >0. The 
calculated <j s h / '(e/167r) is plotted in Fig. 1 as a function 
of a/k F , for three values of V 2 / /3 2 k F =10- 4 , HT 3 , and 
10 _1 . The ratio a/k F is a measure of relative strength 
of the linear to cubic terms in J2J. As expected, the 
<7 s h vanishes for large a. It is important to notice the 
singular behavior at small T of <j s h in the vicinity of 
a/k F =l/y/2 and a/k F =0. The singularities appear be- 
cause at these points the spin-orbit splitting 2hk van- 
ishes for certain k directions. As a result, in such an- 
gular integrals .2^ x — >oo when the elastic scattering rate 
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FIG. 1: Spin-Hall conductivity as a function of a/kF for 
Y 2 /0 1 k%^lQ- i , HT 3 , and HT 1 . 



r^O. It is also interesting to notice that in the range 
0<a/fc F <l/V2, as T^O the spin-Hall conductivity has 
a plateau shape with the universal value of <j s h =3e/87r. 
This plateau and the sharp change of sign at a/k F =l/V2 
can be useful in device applications. 

We would like to elaborate the non-analytic behavior 
of l|lt)fl when both f2 and Q approach zero, a consequence 
of the diffusion denominator in J. When Q^Q first, i|16|) 
gives the DC flow of the spin-Hall current induced by 
the spatially homogeneous electric field. At the opposite 
regime DQ 2 ^>fl, we neglect the fl in the denominator 
and rewrite i|16|) in a coordinate independent form as 



J, = 



iCla. 



H 



DQ 2 



(17) 



Here Ji is the z-polarized spin current flowing along the 
I axis. To arrive at l|17f) we have assumed that DQ 2 is 
much less than the spin relaxation rate. Otherwise, the 
term DQ 2 should be added to (|T3|) . 

i|17|) yields the hydrodynamics of the spin-Hall current 
flow. Since E=-iQV is a longitudinal field, we have 



V • J = ; (V x J) z = 



<J s h dV 
D dt 



(18) 



The first equation indicates that the spin current is con- 
serving. The second equation tells us that in each spatial 
point the flux is perpendicular to the local electric field, 
similar to the spin-Hall effect in a homogeneous field. 
In the field of spherically symmetric potential a circular 
vortex flow of the spin current is thus induced around 
a central charge. The physics of this effect is similar to 
the screening of scalar potential by electric charges. To 
clarify this analogy, let us introduce the conjugate cur- 
rent J x =J y and Jy=-J X i as well as the vortex "charge" 
density p defined by the continuity equation 

eVJ=|. (19) 
We can then rewrite the second equation in l|18|) as 

P = e a -§V, (20) 
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which has the same form as the equation for the electro- 
static screening of the scalar potential V, with ea s njD 
playing the role of the inverse screening length. 

It should be noticed that because of the above men- 
tioned close relationship between the spin-Hall effect and 
the accumulation of in-plane spin polarization, the latter 
will also appear as a screening cloud around the exter- 
nal charge. The in-plane polarization, in its turn, can 
give rise to a z-polarized component via the spin pre- 
cession term of the diffusion equation 30 . This precession 
is proportional to vpQ/T, which is small in the diffu- 
sion approximation and was neglected in (|18|) . Conse- 
quently, the spin-Hall current turns out to be conserved, 
as one can expect in the absence of the relaxation of z- 
polarization. On the other hand, in the near vicinity of 
the vortex core, the precession term becomes more im- 
portant because of the larger gradient of the electric field. 
Hence, the accumulation of the z-polarizcd spin density 
will be expected in the region of the core. The detailed 
analysis of this phenomenon is outside the scope of the 
present paper. It is worthwhile to notice that the core 



has a macroscopic size about Hvf/A, which is of the or- 
der microns. Therefore, the spin accumulation in the 
vortex core can be observed by, for example, the method 
of Faraday rotational. 

In conclusion, within the quasiclassical perturbation 
theory we have shown that, in agreement with existing 
results, impurity scattering reduces the DC spin-Hall cur- 
rent to zero if the SOI is due to the Rashba interaction. 
On the other hand, the spin-Hall current remains finite 
for the Dresselhaus SOI. Nevertheless, this current be- 
comes zero if it is induced by a spatially varying DC elec- 
tric field. The field must be time dependent in order to 
produce a finite effect. In this case the spin-current flow 
in the field of a scalar potential has the form of a vortex. 
The physics of this phenomenon is formally equivalent to 
the screening of external electric potential by electrons. 
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